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Abstract 

Extensive numerical integration results lead us to conjecture that the silver mean, that is, o~Ag 
\ 2 1 .414214 plays a fundamental role in certain geometries (those given by monotone metrics) 
imposable on the i5-dimensional convex set of two-qubit systems. For example, we hypothesize 
that the volume of separable two-qubit states, as measured in terms of (four times) the minimal 
monotone or Bures metric is -jp, and W&Ag in terms of (four times) the Kubo-Mori monotone 
metric. Also, we conjecture, in terms of (four times) the Bures metric, that that part of the 14- 
dimensional boundary of separable states consisting generically of rank-/our 4x4 density matrices 
has volume ( "hyperarea" ) 55 ^ g , and that part composed of T&nk-three density matrices, 3 ^ g , 
so the total boundary hyperarea would be 98 ^ 9 . While the Bures probability of separability (~ 
0.07334) dominates that (« 0.050339) based on the Wigner-Yanase metric (and all other monotone 
metrics) for rank-four states, the Wigner-Yanase (~ 0.18228) strongly dominates the Bures (~ 
0.03982) for the rank-three states. 
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I. INTRODUCTION 



A. Background 

An arbitrary state of two quantum bits (qubits) is describable by a 4 x 4 density matrix 
(-D4) — an Hermitian, nonnegative definite matrix having trace unity. The convex set of all 
such density matrices is i5-dimensional in nature . Endowing this set with the statisti- 
cal distinguishability (SD) metric (identically four times the Bures [minimal monotone] 
metric jjj), we addressed in Q] the question (first essentially raised in the pioneering study 
0| , and investigated further in (f| Q, Qj) of what proportion of the 15-dimensional convex 
set (now a Riemannian manifold) is separable (classically correlated) in nature [jj. This 
pertains to the question of manifest interest "Is the world more classical or more quantum?" 
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nn 

The Peres- Horodecki partial transposition criterion jlOL provides a convenient neces- 



sary and sufficient condition for testing for separability in the cases of qubit-qubit (as well 
as qubit-qutrit) pairs [3]. That is, if one transposes in place the four 2x2 blocks of D4, 
then in the case that the four eigenvalues of the resultant matrix are all nonnegative — or 
more simply, if its determinant is nonnegative Thm. 5] — D4 itself is separable. 

Sommers and Zyczkowski Il4l eq. (4.12)] have recently established (confirming en passant 
certain conjectures of Slater |15j ) that the Bures volume of the (N 2 — l)-dimensional convex 
set of complex density matrices (Djv) of size N is equal to 

2 l-W 2 7r iV 2 /2 



r(N 2 /2) ' 

For the only specific case of interest here, N = 4, this gives us for the total Bures volume 



(1) 



7T 8 



Vg+" = « 5.74538 • 10" 5 . (2) 

Bures 165150720 V ' 

We let the superscript s denote the set of separable and the superscript n the (complemen- 
tary) set of nonseparable 4x4 density matrices. (The comparable volume based on the 
Hilbert- Schmidt metric — which induces the flat, Euclidean geometry into the set of mixed 
quantum states — is = 851 ^ Q5Q0 « 1.12925 ■ 10~ 6 Q eq. (4.5)14 The volume V m?es is 
exactly equal to that of a 15-dimensional haljsphere with radius | 



14j . Now, additionally, 



So, Vgp n is itself exactly equal to one-half the volume ("surface area") of a 15- dimensional 
sphere of radius 1. (The full sphere of total surface area = 2V^ n s ^ s m -^dimensional 
Euclidean space and bounds the unit ball there.) 

One of the objectives in this study will be to highly accurately estimate the included 
volume Vg D . Then, we could, in turn, obtain a good estimate of the SD/Bures probability 
of separability. 
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Also, we could gain evidence as to possible exact values, which on the basis of previous 
lower-dimensional analyses |8[, we have been led to believe is a distinct possibility. 

We had already undertaken this task in [4j (seeking there to exp 
developed Euler angle parameterization of the 4x4 density matrices 
was, however, in retrospect, based on a relatively small number (65 million) of points, gen- 
erated in the underlying quasi-Monte Carlo procedure (scrambled Halton sequences) (cf. 



oft the then just- 
The analysis 



18, 



191]). (Substantial computer assets were required, nonetheless. Numerical integration 
in high- dimensional spaces is a particularly challenging computational task.) One of the 
classical "low-discrepancy" sequences is the van der Corput sequence in base b, where b is 
any integer greater than one. The uniformity of the van der Corput numbers can be further 
improved by permuting/scrambling the coefficients in the digit expansion of N in base b. 



The scrambled Halton sequence in iV-dimensions — which we employed in j4j and in our 



auxiliary analyses below (sec. IIV|) — is constructed using t 
numbers for 6's ranging over the first N prime numbers [19 



re so-scrambled van der Corput 
p. 53]. 



To facilitate comparisons with the results of Sommers and Zyczkowski [14], which were 
reported subsequent to our analysis in 0, we need to both divide the estimates given in jj] 
by 4! = 24 to take into account the strict ordering of the four eigenvalues of D4 employed by 
Sommers and Zyczkowski jl^, eq. (3.23)], as well as to multiply them by 8, since we (due to 
a confusion of scaling constants) only, in effect, used a factor of 2 12 in |4j, eq. (5)-(7)] rather 
than one of 2 15 , as indicated above in (J3J is required. These two independent adjustments 



s+n 



together amount to a multiplication by 8/24 = 1/3. This means that the estimate of V SD 
(the true value of which, as given above, is known to be ~ 1.882645) from the quasi-Monte 
Carlo analysis in jj], should be taken to be 1.88284 = 5.64851/3; the estimate of Vg D from 
P] should, similarly, be considered to be 0.138767 = .416302/3; and of P§ D (for which no 
adjustment is needed, being a ratio), 0.0737012. 



We had been led in J^jj] — if only for numerical rather than any clear conceptual reasons 
- to formulate a conjecture that (after adjustment by the indicated factor of |) can be 
expressed here as 

VS n = = 0.138729, (5) 

SD 6g30 , V ) 

as well as that 

Psd = P Lres = 7^™ 0.0736881 (6) 
(suggesting that the [quantum] "world" — even in the case of only two qubits — is consider- 
ably "more quantum than classical"). We now must view (JSJ) and (jUJ) as but approximations 
to the revised conjectures ()15|) and ()16|) below, obtained on the basis of such much larger 
quasi-Monte Carlo calculations. 



B. Monotone Metrics and Quasi-Monte Carlo Procedures 



The Bures metric plays the role of the minimal monotone metric. The monotone metrics 



comprise an infinite (nondenumerable) class 
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221 ]. generalizing the (classically unique) 



Fisher information metric |2 
member of this class 



he Bures metric has certainly been the most widely-studied 



For the infinitesimal distance element between 



two states -D4 and -D4 + SD4, we have 



(ds B i 



1 



(?) 



where D4 is diagonal in the orthonormal basis {\j)} with eigenvalues (A?}. 

Two other prominent members are the maximal monotone metric 1 281 and the Kubo-Mori 
3iJ,|3Jj (also termed Bogoliubov-Kubo-Mori and Chentsov [321 ]) monotone metric. 
The Kubo-Mori metric (or canonical correlation stemming from differentiation of the relative 
entropy) is, up to a scale factor, the unique monotone Riemannian metric with respect to 
which the exponential and mixture connections are dual [2^], and as such, certainly merits 
further attention. 

In this study, we will utilize additional computer power recently available to us, together 
with another advanced quasi-Monte Carlo procedure (scrambled Faure-Tezuka sequences 
[33^ - - the use of which was recommended to us by G. Okten, who provided a corre- 
sponding MATHEMATICA code). Faure and Tezuka were guided "by the construction 

C (i) = ^(i)p(i-i) and by 

some possible extensions of the generator formal series in the 



framework of Neiderreiter" . (A^ is an arbitrary nonsingular lower triangular [NLT] matrix, 
P is the Pascal matrix j^J and is a generator matrix of a sequence X). Their idea 
was to multiply from the right by nonsingular upper triangular (NUT) random matrices 
and get the new generator matrices = p( l-1 )[/W for (0, s)-sequences [33] . "Faure-Tezuka 
scrambling scrambles the digits of % before multiplying by the generator matrices . . . The 
effect of the Faure-Tezuka-scrambling can be thought of as reordering the original sequence, 
rather than permuting its digits like the Owen scrambling . . . Scrambled sequences often have 



smaller discrepancies than t 
facilitate error estimation" 



reir nonscrambled counterparts. Moreover, random scramblings 



m p. io7]. 



The Faure-Tezuka procedure appears to us to be exceptionally successful in generating 
a highly uniform (low discrepancy j^]) distribution of points over the hypercube - - as 
judged by its yielding an estimate of 1.88264 for V^ n ~ 1.882645. However, at this stage, 
the procedure does have the arguable shortcoming that it does not readily lend itself to 
the use of "error bars" for the estimates it produces, as quite naturally do (the generally 
considerably less efficient) Monte Carlo methods (which, of course, distribute points on the 
basis of pseudorandom, rather than deterministic, methods.) 

"It is easier to estimate the error of Monte Carlo methods because one can perform a 
number of replications and compute the variance. Clever randomizations of quasi-Monte 
Carlo methods combine higher accuracy with practical error estimates" j^j], p. 95]. G. 
Okten is presently developing a MATHEMATICA version of the scrambled Faure-Tezuka 
sequence in which there will be a random generating matrix for each dimension — rather 
than one for all [fifteen] dimensions — which will then be susceptible to statistical testing 

0. 

C. Morozova-Chentsov Functions 



To study such monotone metrics other t 



ran the SD/Bures one, we will utilize a certain 



ansatz (cf. (3|). Contained in the formula 3, eq. (3.18)] of Sommers and Zyczkowski for 
the "Bures volume of the set of mixed quantum states" is the subexpression (following their 
notation), 

n tj1...N (Pv ~ Pm) /o\ 
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where p M , p u (//, v = 1, . . . , N) denote the eigenvalues of an N x N density matrix (-Djv)- 
The term (jHJ) can equivalently be rewritten using the "Morozova-Chentsov" function for the 
Bures metric eq. (2.18)], 

2 

C B ures(Pfi, Pu) = ; , (9) 

Pu + P» 



as 



Qn = (p v - p^) 2 c Bures (p^ p v )/2. (10) 

A Morozova-Chentsov function is a positive continuous function c(A,/i) that is sym- 
metric in its two variables and for which c(A, A) = CA -1 , for some constant C, and 

n 

c(t\,tp) = t- l c(\p) [21, Thm. 1.1]. There exist one-to-one correspondences between 



Morozova-Chentsov functions, monotone metrics and operator means. |2JJ, Cor, 6]. "Oper- 
ator means are binary operations on positive operators which fulfill the main requirements 
of monotonicity and the transformer inequality" j^ . 

The ansatz we employ is that the replacement of cs ures (p M , p u ) in the formulas for the 
Bures volume element by the particular Morozova-Chentsov function corresponding to a 
given monotone metric (g) will yield the volume element corresponding to that particular 
g. We have been readily able to validate this for a number of instances in the case of the 
two-level quantum systems [N = 2L using the general formula for the monotone metrics over 



such systems of Petz and Sudar |2Qj, eq. (3.17)]. One can argue that the joint distribution of 
the eigenvalues of D N is the product of Qn — pertaining to the off-diagonal elements of the 
density matrix — and an additional factor — pertaining to the diagonal elements. Now, 
H jsj is equal to the reciprocal of the square root of the determinant of the density matrix 
for all [Fisher-adjusted] monotone metrics — so we need not be concerned with its variation 
across metrics in this study — and simply unity in the case of the [flat] Hilbert-Schmidt 
metric (cf. 

D. Outline of the study 

In addition to studying the SD/Bures metric, we ask analogous questions in relation to 
a number of other monotone metrics of interest. We study two of these metrics, in addition 
to the SD metric, in our "main analysis" (sec. IIII|) and two more in our "auxiliary analysis" 
(sec. IIV|) . which is based on the same scrambled Halton procedure employed in ^] — but 



with more than five times the number of points generated there, but also many fewer points 
than in the primary (main) analysis here. (In hindsight, we might have better consolidated 
the several monotone metrics into a single investigation, from the very outset, but our 
initial/tentative/exploratory analyses grew, and we were highly reluctant to discard several 
weeks worth of demanding and apparently revealing computations. Also, we had been using 
two different sets of processors [Macs and Suns] for our computations and for a number 
of reasons — too involved and idiosyncratic to make the subject of discussion here — it 
proved convenient to conduct two distinct analyses.) Also, we include analyses in sec. IS 
pertaining to the maximal monotone metric, and a number of metrics interpolated between 
the minimal and maximal ones. (The "average" monotone metric — studied in our main 
analysis (sec. I111J) — is obtained by such an interpolation.) In sec. IV II we apply Monte-Carlo 
methods in a limited study of the questions raised before. In sec. I VI II we undertake studies 
concerned with the values of volumes ( "surface areas" ) of the 14- dimensional boundary of the 
15-dimensional convex set of two-qubit states, as measured in terms of the various monotone 
metrics under investigation here. 

To begin with (sec.|n}, we will seek to determine V^Q- A wiggly line over the acronym for 
a metric will denote that we have ab initio multiplied that metric by 4, in order to facilitate 
comparisons with results presented in terms of the SD, rather than the Bures metric, which 
is one-fourth of the SD metric. (This, perhaps fortuitously, gives us a quite appealing scale 
of numerical results.) The probabilities themselves — being computed as ratios — are, of 
course, invariant under such a scaling, so the "wiggle" is irrelevant for them. 



II. PRELIMINARY ANALYSIS OF THE KUBO-MORI METRIC 



The Morozova-Chentsov function for the Kubo-Mori metric is [14, eq. (2.18)] 

logp„ -logp M . . 

ckm[Ph,Pu) = -• (11) 

Pv PjX 

To proceed in the study of the KM metric, we first wrote a MATHEMATICA program, 
using the numerical inte grat ion command, that succeeded to a high degree of accuracy in 
reproducing the formula jlj, eq. (4.11)], 

2 A^r(A72) 

07V " 7HW(i)...r(;v + i) 1 j 
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TABLE I: Estimates based on (four times) the Bures, "average" and Kubo-Mori metrics, us- 
ing a scrambled Faure-Tezuka sequence composed of two billion points distributed over the 15- 
dimensional unit hypercube, for quasi-Monte Carlo numerical integration. The results based on 
the first one billion points are given in parentheses. 



metric 


ys+n 

•metric 


Metric 


P s . . =v s iv s+n 

metric 1 


Bures 


1.88264 (1.88264) 


0.137884 (0.137817) 


0.0732398 (0.0732042) 


Average 


28.0801 (28.0803) 


1.33504 (1.33436) 


0.0475438 (0.0475194) 


Kubo-Mori 


120.504 (120.531) 


4.1412 (4.14123) 


0.0343654 (0.0343583) 



for the Hall/Bures normalization constants [bol for various N . (These constants form 
one of the two factors — along with the volume of the flag manifold [1J, eqs. (3.22), (3.23)] 
- in determining the total Bures volume.) Then, in the MATHEMATICA program, we 
replaced the Morozova-Chentsov function Q for the Bures metric in the product formula 
(|10J) by the one ()11|) for the Kubo-Mori function. For the cases A^ = 3, 4 we found that the 
new numerical results were to several decimal places of accuracy (and in the case N = 2, 
exactly) equal to 2 N ( N ~ 1 ^ 2 times the comparable result for the Bures metric, given by (fT2*|) . 
This immediately implies that the KM volumes of mixed states are also 2 N ^ N ~ 1 ^ 2 times 
the corresponding Bures volumes (and the same for the KM and SD volumes), since the 
remaining factors involved, that is, the volumes of the flag manifolds are common to both 
the Bures and KM cases (as well as to all the monotone metrics). Thus, we arrive at our 
first conjecture (cf. jl4j|). 

4vr 8 



ys+n _ aAys+n 
KM ~ ^ V SD 



315 



120.489, 



(13) 



for which we will obtain some further support in our main numerical analysis (sec. |TTT 
yielding Table |U 



III. MAIN ANALYSIS 



Associated with the minimal (Bures) monotone metric is the operator monotone function, 
fn,, rps (t) = (1 + t)/2, and with the maximal monotone metric, the operator monotone 
Lion, fl (t) I mi + t) Q e, (2,7)]. The me ra S e of the8 e two taction, that is, 



s 



faverage(t) = (1 + 6t + t 2 ) / (4 + it), is also necessarily operator monotone [2l|, eq. (20)] and 
thus yields a monotone metric (apparently previously uninvestigated). Again employing our 
basic ansatz, we used the associated Morozova-Chentsov function — given by the general 
formula [20j, p. 2667], c(x,y) = l/yf(x/y) — 

_ 4(p M + p v ) 

Cavera9e[P " ,Pu) ~ Pl + Gp»p„ + Pl- 1 J 

For our main quasi- Monte Carlo analysis, we (simultaneously) numerically integrated the 
SD, KM and avg volume elements over a fifteen-dimensional hypercube using two billion 
points for evaluation, with the points forming a scrambled Faure- Tezuka sequence j^J • (As 
in 4j , the fifteen original variables — twelve Euler angles and three angles for the eigenvalues 
|l7l eq. (38)] — parameterizing the 4x4 density matrices were first linearly transformed so 
as to all lie in the range [0,1].) This "low-discrepancy" sequence is designed to give a close- 
to- uniform coverage of points over the hypercube, and accordingly yield relatively accurate 
numerical integration results. 

The results of Table U suggest to us, now, rejecting the previous conjecture © — based 
on a much smaller number (65 million) of data points than the two billion here — and 
replacing it by (perhaps the more "elegant") 

V§ D = ^f = « 0.138071, (15) 



where o~Ag denotes the "silver mean" [38]. (As we proceed from one billion to two billion 
points, some apparent convergence — 0.137817 to 0.137884 — of the numerical estimate to 
the conjecture (fT5|) is observed. It is interesting to note the occurrence of the first three 
positive integers in (|TK|) — a property which obviously the much-studied golden mean, 
lacks.) By implication then, the conjecture (jOJ) is replaced by 



P S so,Bur es = ^ = 1 ~^ L - 0.0733380. (10) 

V SD n 



In addition to simply our numerical results, we were also encouraged to advance this 

n 

conjecture (|15|) on the basis of certain earlier results. In |8[, a number of quite surprisingly 
simple exact results were obtained using symbolic integration, for certain specialized [low- 
dimensional] two-qubit scenarios. This had led us to first investigate in ^| the possibility of 
an exact probability of separability also in the full 15-dimensional setting. (Unfortunately, 
as far as we can perceive, the full 15-dimensional problem is not at all amenable — due to 



its complexity — to the use of the currently available symbolic integration programs and, it 
would appear, possibly for the foreseeable future.) 

In particular in |sjj], a Bures probability of separability equal to a^ g had been obtained for 
both the q — 1 and q — | states inferred using the principle of maximum nonadditive 
[Tsallis] entropy — and also for an additional /otw-dimensional scenario (s| sec. II.B.l]. (We 
have recently reanalyzed this last scenario, but with the maximal monotone metric, and also 
found a probability of separability equal to a^ g . The value a^g also arises as the amount 
by which Bell's inequality is violated 4(J eq, (8)].) Christos and Gherghetta 38] took the 
silver mean, as in our study here, to be the positive solution of the equation x + 2 = - 



since having a "mean" value less than 1 was use 



scaling functions — while others (perhaps more) 



41 



ul in their investigation of trajectory 



Q 



chap. 22] have defined it as 



the positive solution of x — 2 = 1, that is, v^2 + 1, the reciprocal of our o~Ag- (The square root 
of two minus one is also apparently a form of "Pisot number" 44j. Similar definitional, but 
perhaps not highly signficant, ambi guit ies occur in the (more widespread) usage of the term 
"golden mean", that is (a/5 ± l)/2 [4.51 • ("The characteristic sequence of (V5-l)/2 (resp., 
y/2 — 1) is called the golden mean sequence (resp., Pell sequence)" This line of analysis 
- concerned with the alignment of two words over an alphabet — originated, apparently, 
from a 1963 unpublished talk of the prominent [Nobelist] physicist, D. R. Hofstadter j3|.) 
In j47|, demonstrating a conjecture of Gromov, the minimal volume of R 2 (the infinite 
Euclidean plane) was shown to be Jj-. (An exposition of this result is given in 4^1-) ^ n 3 
the value of was obtained for a certain supremum of volumes. 

Further conjectures that V& g = « 1.33469 and V S R ~ M = 10a Ag » 4.14214 seem worth 
investigating, based on the results in Table |l] (Our estimate of the ratio -r^ 1 from Table H] 

V SD 

is 30.0339.) So, we have an implied conjecture that 

Pkm = V M = - -0343776. (17) 

KM Z7F 

It would then follow that p / kM =§ = .46875. 

SD j Bures 

The convergence to the known value of V^ n in Table H] seems more pronounced than any 
presumptive convergence to the conjectured values of the separable volumes alone, but the 
latter are based on considerably smaller samples (roughly, one-quarter the number) of points 
than the former (for which, of course, all the two billion systematically generated points are 
used). Clearly, our conjecture (fTHj) can be reexpressed as Vg D = V^ n — V§ D = — -da. m 
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FIG. 1: Deviations of the estimated values of V^ n from the known value — as shown by Sommers 
and Zyczkowski — of ~ 1.882645, as the number of points in the scrambled Faure-Tezuka 
sequence increases from 1 million to 2,000 million 




FIG. 2: Deviations of the estimated values of V^~™ from the conjectured value of ~ 120.489, 
as the number of points in the scrambled Faure-Tezuka sequence increases from 1 million to 2,000 
million 

1.74457. Our sample estimate for Vg D is, then, 1.74475.) 

In Figs. WS^ we show the deviations from our conjectured and known values of the es- 
timates provided by the Faure-Tezuka sequence as the number of points in the sequence 
increases from one million to two thousand million (i. e. two billion). 

In Fig. El additionally, we show together the relative deviations of V^ n and of from 
the known and conjectured values. In other words, we divide the estimated values by the 
known/conjectured values and subtract 1. The SD curve is extraordinarily better behaved 
("hugging" the x-axis) than is the KM curve. Perhaps this difference is attributable to the 
"simpler" (more numerically stable?) nature of the Morozova-Chentsov function in the SD 
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FIG. 3: Deviations of the estimated values of Vg D from the conjectured value of — ^ ~ 0.138729, 
as the number of points in the scrambled Faure-Tezuka sequence increases from 1 million to 2,000 
million 
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FIG. 4: Deviations of the estimated values of V^ g from the conjectured value of — g- a ~ 1.33469, 
as the number of points in the scrambled Faure-Tezuka sequence increases from 1 million to 2,000 
million 

case © than in the KM case (ITU . 

Further plotting of our various results yielded one of particular interest. In Fig. we 
show the estimates of V^™ 9 — Vg D . Of course, this figure (the scale of which was internally 
chosen by MATHEMATICA, based on the data, and not exogeneously imposed) strongly 
suggests that V a " s — Vg D = 25. Now, we found that if we posit 

2 C )7T 8 

V$T = « 28.0792, (18) 

(with 8448 = 2 8 • 3 • 11), it would follow that 

244Q7T 8 26 

~ V s n D = - -a Ag » 24.99996094, (19) 
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FIG. 5: Deviations of the estimated values of from the conjectured value of Wo~A g ~ 4.14214, 
as the number of points in the scrambled Faure-Tezuka sequence increases from 1 million to 2,000 
million 
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FIG. 6: Joint plot of relative deviations of estimates of V^ n and from their known and 

conjectured values of and |^-. The more rugged curve corresponds to V^Q- 

(with 887040 = 2 8 • 3 2 • 5 • 7 • 11), being strikingly close to the indicated value of 25. 

In Fig. IHlwe plot the deviations of the estimates of V^ n from its conjectured value (fl~8|) . 



IV. AUXILIARY ANALYSIS 



In an independent set of computations (Table ITT|). employing 415 million points of a scram- 
bled Halton sequence (as opposed to the scrambled Faure-Tezuka sequence used in sec. IHIj) . 
we sought to obtain estimates of the probability of separability of two arbitrarily coupled 
qubits based on three monotone metrics of interest. The specific method employed was 
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FIG. 7: Estimates of the difference, V£~ g — Vg D , as the number of points in the scrambled Faure- 
Tezuka sequence increases from 1 million to 2,000 million. Note the strong suggestion that the true 
value is 25 (or close thereto) 
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FIG. 8: Deviations of the estimated values of V^~t n from the conjectured value of §§fg , as the 
number of points in the scrambled Faure-Tezuka sequence increases from 1 million to 2,000 million 

that of scrambled Halton sequences [3]. (While there are different scrambled Faure-Tezuka 
sequences depending upon the particular random generating matrices used, the scrambled 
Halton sequence is unique in nature.) 

These correspond to the operator monotone functions, 



fGKs{t)=t t l^/t 



fwY(t) = ^(Vt+l) 2 ; f KM (t) 



t - 1 
logt 



(20) 



The subscript GKS denotes the Grosse-Krattenthaler-Slater ( "quasi-Bures" ) metric (which 
yields the common asymptotic minimax and maximin redundancies for universal quantum 



coding 
sec. 4] jo 



501 sec. IV.B1 



5l(), the subscript WY, the Wigner-Yanase information metric 



52, 



, and the subscript KM, the Kubo-Mori metric already studied in sees. II and 
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TABLE II: Estimates based on (four times) the Grosse-Krattenthaler-Slater, Wigner-Yanase and 
Kubo-Mori monotone metrics, using a scrambled Halton sequence consisting of 415 million points. 



metric V s+ ? . V s -. P* tric = V s /V s 

metric metric men ic ' 

GKS 5.4237 0.330827 .0609965 

WY 14.5129 0.730567 .0503391 

KM 120.504 4.1791 .0346801 



III. (We had, in fact, intended to study the "Noninformative" monotone metric |55J] here 
instead of the KM metric, but there was a programming oversight that was only uncovered 
at the end of the computations.) It appears conjecturable that, in terms of the separable 
states, Vqk S = ~ 0.331371. (The evidence is somewhat of a weaker nature that 

V^y = ~I S ' ~ 0.724874.) Also, in terms of the combined separable and nonseparable 
states, it seems possible that V£g s = ^ » 5.42202, with 1750 = 2 ■ 5 3 ■ 7. If so, we would 
have P^ KS = ±^d£ « 0.0611158. 

We had hoped to further extend the scrambled Halton sequence used here, but doing so 
has so far proved problematical, in terms of available computer resources. 



V. MAXIMAL MONOTONE METRIC 



As to the maximal monotone metric, numerical, together with some analytical evidence, 
strongly indicate that is infinite (unbounded) (as well as V^ ax ). The supporting ana- 
lytical evidence consists in the fact that for the three-dimensional convex set of 2 x 2 density 
matrices, parameterized by spherical coordinates [r, 9, <fi\ in the "Bloch ball", the volume 
element of the maximal monotone metric is r 2 sin 9(1 — r 2 )~ 3 / 2 , the integral of which di- 
verges over the ball. Contrastingly, the volume element of the minimal monotone metric is 
r 2 sin#(l — r 2 ) -1 / 2 , the integral over the ball of which is finite, namely tc 2 . For s > 1 the 
integral of r 2 sin 9(1 — r 2 )~ s diverges, so the divergence associated with the monotone metric 
itself is not simply marginal or "borderline" in character. 

To gain further evidence in these regards, one can engage in numerical estimation for the 
one-parameter family of interpolating metrics given by the operator monotone functions 

f a (t) = (1 - a)f max (t) + af Bures (t), (21) 
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TABLE III: Estimates based on the first ninety-six million points of a scrambled Faure-Tezuka 
sequence of (four times) a number of metrics obtained by interpolating between the maximal 
(a = 0) and minimal (a = 1) monotone metrics 

a 1 KT 1 10" 2 10" 3 10" 4 10" 5 

V* +n 1.88258 7951.27 

V a s 0.13786 148.569 

P a s 0.073229 0.01868 



9.3254 • 10 b 6.0345 • 10 9 3.049 • 10 12 1.2825 • 10 i5 8.0858-10* 
63659. 2.0972 • 10 7 6.502 • 10 9 2.2084 • 10 12 5.229 • 10 36 
0.006827 0.003475 0.0021325 0.0017219 0.00064669 



for which the Morozova- Client sov functions are of the form 

c a {p^Pu) = —, \ 2 , , • 22 

a(p M - p v y + Ap^p u 

Then, one could plot the results as a function of the parameter a and study the limit a — > 0. 
(Of course, for a = h, one would recover the "average" monotone metric, studied in sec. 11111 ) 
A preliminary investigation along these lines is reported in Table IIHI Based on the first 
ninety-six million points of a scrambled Faure-Tezuka sequence, we obtain estimates of V^ +n , 
and P^ = J s a +n for (four times) the metrics interpolating between the maximal (a = 0) 
and minimal (a = 1) monotone metrics for several values of a, increasingly close to a = 0. 
So, P^ax wou ld seem quite close to being 0. (However, some clearly numerically anomalous 
behavior occurred in passing from ninety-six million points to ninety-seven million points. 
The estimates of V*£, V^ ax and P s max jumped to 8.27999 ■ 10 40 , 7.48026 ■ 10 40 and 0.903414, 
respectively.) 

It would be interesting to formally test the hypothesis that P^ ax = 0. (More specifically, 
we might ask the question if the limit of P£ as a — > is 0.) If it can, in fact, be established 
that P^ax i s zero ) this might serve as something in the nature of a "counterexample" to 
the proposition (a matter of considerable interest in the theoretical analysis of quantum 
computation) that for bipartite quantum systems of finite dimension, there is a separable 



neighborhood of the fully mixed state of finite volume [56|, |57|, |58j, |59j . (These conclusions 
were obtained with the use of either the trace or Hilbert-Schmidt metric — the first of which 
is monotone, but not Riemannian, while the second is Riemannian, but not monotone |x3| .) 

We have conducted a test along these lines. Using a simple Monte- Carlo (rather than 
quasi- Monte Carlo) scheme, we generated ten sets of ten million points randomly distributed 
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over the 15- dimensional hypercube. For each of the ten sets, we obtained estimates of V^+™, 
V^ax and hence P max - Based on the one hundred million points the (mean) estimate of P max 
was /i = 1.77038 • 10~ 7 and the standard deviation across the ten samples, rj = 3.692 • 10~ 7 . 
So, the value lies less than one-half (that is, 0.479510) standard deviations from \i. For a 
student t-distribution with 9 = 10 — 1 degrees of freedom, forty percent of the probability 
lies outside 0.261 standard deviations from the mean and twenty- five percent outside 0.703 
standard deviations. So there is little evidence here for rejecting a hypothesis that P max 
equals 0. For an independent analysis based on ten sets of four million points, the estimates 
were roughly comparable, i. e., /i = 2.4196 • 10~ 7 , 77 = 5.09683 • 10~ 7 . Also, for ten sets of 
five million points, but setting the interpolation parameter a not to but to .05, there were 
obtained fi = .00438593 and 77 = .000229437, with H = 19.1611. So, here one can decisively 
reject a hypothesis that the probability of separability for a = .05 is 0. 



VI. FURTHER MONTE-CARLO ANALYSES 



We also undertook a series of Monte- Carlo analyses, incorporating together the GKS, Bu- 
res, average, Kubo-Mori, Wigner-Yanase, maximal and Noninformative (NI) 5^| monotone 
metrics. (The operator monotone function f{t) associated with the NI metric is n_^7og 2 (t) ') 
We subdivide the unit hypercube into 3 15 = 14, 348, 907 subhypercubes, pick a random 



point in each one of these, and then repeat the procedure... We are now able to report in 
Table IIVI the results of fifteen iterations of this process. The central limit theorem tells us 
that for a large enough sample size, the distribution of the sample mean will approach a 
normal/Gaussian distribution. This is true for a sample of independent random variables 
from any population distribution, so long as the population has a finite standard deviation. 
The population standard deviation is equal to the standard deviation of the mean times 
the square root of the sample size N, which in our case is 15 • 3 15 . If one were to use two 
standard deviations as a rejection criterion, then the only one of our conjectures that would 
be rejected would be that for V^„. (However, the standard deviations in the separable 
cases would be approximately four times as large if we only used the number of points cor- 
responding to separable states, rather than to all states, as we have done here. This would 
lead, then, to none of our conjectures being rejected.) The estimate of P max , obtained in 
the same Monte-Carlo procedure, was jf|f§^ = 2.24073 • lO" 10 . 
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TABLE IV: Monte-Carlo analysis based on 15 • 3 = 215, 233, 605 density matrices. The second 
set of columns correspond to V s - . and the third, to V s+ ~ n . . Estimates (est.) of the volumes, 

metric metric v ' 



standard deviations (s.d.) of these estimates and the number of standard deviations they are away 
from their conjectured (cj.) or known values, as given in Table IVfl are presented 



metric 


est. 


s. d. 


cj.—est. 
s.d. 


est. 


s. d. 


cj.—est. 
s.d. 


Bures 


0.13800 


0.00021 


0.3354 


1.88295 


0.00102 


-0.2985 


GKS 


0.32990 


0.00057 


2.5696 


5.4232 


0.00315 


-0.3867 


WY 


0.72811 


0.00152 




14.5084 


0.00943 




Avg 


1.3363 


0.00287 


-0.5821 


28.0781 


0.01769 


0.0618 


KM 


4.1574 


0.02242 


-0.6792 


120.256 


0.17489 


1.3363 


NI 


848.05 


28.997 




48668.2 


421.712 





VII. ^-DIMENSIONAL BOUNDARIES 



A. Initial analyses 

In the analyses above, we have been concerned with the volume of the 15- dimensional 
convex set of 4 x 4 density matrices, as measured in terms of a number of monotone metrics. 
We have modified the computer programs involved, so that they would provide estimates of 
the volume ( "hyper area" ) of the boundary of this set. 

Our numerical integrations were conducted over a /emrteen-dimensional hypercube, now 
allowing one of the original fifteen variables (specifically, the hyperspherical angle designated 
#3 in (j, eq. (2)]) to be determined not by the quasi- Monte Carlo procedure, but by the 
requirement that the determinant of the partial transpose equal zero. This considerably 
increases the computational effort per point generated. 

Our early estimates, in this regard, were: 0.587532 (SD), 6.25466 (avg) and 19.8296 
(KM), all three based on the first 2,600,000 points of a scrambled Faure-Tezuka sequence; 
1.47928 (GKS), 3.37384 (WY) and 19.9277 (KM), all three based on the first 1,500,000 
points of a scrambled Halton sequence; and 0.588816 (a=l), 837.072 (a = 10" 1 ), 414676. 
(a = TO" 2 ), 1.57088 • 10 8 (a = 10~ 3 ), 4.87246 • 10 10 (a = KT 4 ), 1.30774 • 10 13 (a = 10~ 5 ), 
8.57852 ■ 10 27 (a = 0) (all seven based on the first 1,100,000 points of a scrambled Faure 
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TABLE V: Estimates of 14-dimensional boundary volumes based on (four times) the Bures, Grosse- 
Krattenthaler-Slater, Wigner-Yanase, Average, Kubo-Mori and Noninformative metrics, using a 
scrambled Faure-Tezuka sequence composed of thirty-five million points. 



metric B s B s+n 

Bures 0.456593 11.4443 0.584072 

GKS 1.45204 

WY 3203.81 17576.3 6606.58 

Avg 6.60067 246.716 6.20592 

KM 19.6215 

NI 4333.36 



sequence). 

Let us note that in terms of the Bures metric — identically one-fourth of the statistical 
distinguishability (SD) metric — the pure state [rank 1] boundary of the 4x4 density 

3 

matrices, both separable and nonseparable, is known to have volume ^- « 5.16771. (This is 
equal to the volume of a 6-dimensional ball of radius 1 and to the volume of a 3-dimensional 
complex projective space |l3 . sec. IV. C].) The 14-dimensional [(3(3] submanifold of 4 x 4 
density matrices of rank 3 has Bures volume 432 ^ 7 32Q ~ 6.98444 • 10~ 4 . Multiplying by 2 14 , 
we obtain the SD counterpart to this of 

51 2ir 7 

B s + n = « 11.4433. (23) 

SD 135135 V ; 

This is twice (a "double-covering") the (n — l)-content (surface area or hyperarea) of the 
unit sphere in n = 15 dimensions. 



B. Further analysis 

Subsequently, we joined all the monotone metrics of interest into a single joint analysis, 
using an independent Faure-Tezuka sequence of points in the 14-dimensional hypercube. 
Up to this point in time, we have generated thirty-five million points (Table IV)). For each 
of these points, we sought values of the fifteenth coordinate — #3 — for which the partial 
transpose of the corresponding 4x4 density matrix had zero determinant. For 24,038,658 of 
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FIG. 9: Deviations of the estimated value of B^jJ 1 from the known value of §g^g ~ 11.4433, as 
the number of points in a certain scrambled Faure-Tezuka sequence increases from one hundred 
thousand to thirty-five million 
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FIG. 10: Deviations of the cumulative estimates of B S SD from the conjectured value of 3 ^ g ~ 
0.456697. The number of points are recorded in steps of 100,000. 

the points at least one feasible value of 9% was found. The even-numbered solutions strongly 
dominated the odd-numbered solutions. (This may pertain to the fact that in this series of 
analyses, we had used — adjusting accordingly — the range [0, 7r], rather than [0, |], as in 
our other analyses, for 9%, 9<i and #3.) There were 74 points with one solution, 30 with three, 
and 11 with five, while there were 2,553,168 with two, 21,312,933 with four, 172,429 with 
six and 3 with eight. 

In Fig. we show the cumulative approximations (in steps of one hundred thousand 
points) to the known value (|2^j) of B%t?. We conjecture (Fig. [TUJ) that the component of 
B s sd consisting of separable states 6l(, that is Bg D , has the value -^p- ~ 0.456697. The 
concomitant estimate of the SD/Bures probability of separability of such rank-three states 
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FIG. 11: Deviations from the conjectured value of (3sd = —^f^ ~ 0.584147 of the cumulative 
estimates of the SD volume of that part of the 14-dimensional boundary of the separable states 
consisting of nondegenerate states. The number of points are recorded in steps of 100,000. 

would then be 

^—=^=?S^ - 0.0398167. (24, 

This is considerably less than the general probability of separability [of, generically, rank- 
four states], conjectured in formula (fboj) to be 0.0733389. The ratio of these two probabilities 
is ^gf « 0.542194. 

In Fig. El we show the cumulative approximations to a conjectured value of (3sd — 
55 ^g 3 ~ 0.584147 for the 14-dimensional SD boundary of separable two-qubit states com- 
posed of generically rank-four states. (The test for membership in this class is that the deter- 
minant of the partial transpose of the corresponding 4x4 density matrix be zero.) Thus, we 
have an implied conjecture that the 14-dimensional boundary of separable 4x4 density matri- 
ces has total SD-volume of (3sd + B s sd = 98 J g Ag « 1.04084. The fit of our cumulative estimates 
to this conjecture is shown in Fig. E3 Numerical evidence (Fig. H*3| also possibly suggests 
that B s a - g = « 6.60153; that (Fig. ED B s Wy = 7735^ « 3203.94; that (Fig. ED 

= 2 9 • 3B S S + D n « 17576.9; that (Fig. ED P K ~ M = « 19.6274; and that (Fig. ED, 

@gks = 27 77° 9 ~ 1.45244. So, surprisingly, the probability of separability of a rank-three 
state appears to be much higher, that is, ^^SiIa^ 9 = 357 17(TAb ~ 0.182281 using 
the Wigner-Yanase metric (but not the average metric, for which we have a sample estimate 
of 0.0267541 and a conjecture of 10 ™g^ j49 ~ 0.0267572) than with the Bures or SD metric, 
in strong contrast to the rank- four case examined earlier (sees. IIIII and lTVj) . 

As is apparent, there are no estimates reported in Table IVl for B s and B s+n for the GKS, 
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FIG. 12: Deviations from the conjectured value of B S SD + (3 S d = ~ 1.04084 of the SD- 

volume of the total boundary of separable states — composed of nondegenerate (rank-four) and 
degenerate (rank-three) 4x4 density matrices. The number of points of the Faure-Tezuka sequence 
are recorded in steps of 100,000. 
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FIG. 13: Deviations from the conjectured value of — ^A^- 6.60153 of the cumulative estimates of 

DS 

avg- 

KM and NI metrics. In retrospect, we could have included the GKS metric (taking the limit 
of the corresponding volume element as one of the eigenvalues of the associated 4x4 density 
matrix approaches zero), but the volume elements for the other two were found to diverge 
on the rank-three states, so they could not have been included. 



C. Levy-Gromov Isoperimetric Inequality 

The scalar curvature of the Bures metric for the 4x4 density matrices is bounded be 
by 570 0, 

Cor. 3]. However, application of the Levy-Gromov Isoperimetric Inequality 
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FIG. 14: Deviations from the conjectured value of 7735aAg ~ 3203.94 of the cumulative estimates 
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FIG. 15: Deviations from the conjectured value of B^±™ of 2Q H^f 
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App. C] requires a lower bound of N on Ricci(Y,Y), where Ricci is the Ricci tensor and 
Y runs over all unit tangent vectors, for closed (N + 1) -dimensional manifolds, We did not 
immediately know if this condition is satisfied or not (given that one can not apparently 
"control" the Ricci curvature in terms of the [bounded] scalar curvature) , but we have found 
that the inequality is violated, and that the condition is not satisfied in the case before us. 

To reach this conclusion, we first took the parameter a (strictly following the notation 
in j^l) to be PsD/Bures = 77§& ~ 0.0736881, according to our conjecture © above. Then, 

V SD 



the function s(a) is the 14- dimensional volume of the boundary sphere dB a , where the 
volume of the ball B a itself is equal to a vol(S 15 ), and 5* 15 is the standard 15-dimensional 
sphere. Further, the function Jsis(a) is the ratio of s(a) ~ 0.499459 to vol(S 



15 a 



2027025 



0.381443. The Levy-Gromov Inequality then asserts that Is\5(a), which here equals 1.30939, 
must be less than a certain ratio, which in our case would — according to our conjectures 
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FIG. 16: Deviations from the conjectured value of (3 K ~, 
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FIG. 17: Deviations from the conjectured value of (3 



GKS 



77 



1.45244 



and known values — be 



B 



Bures 



+ P Bures 2 (Bg D + (3, 



SDj 



Trs+n 

Bures 



z V SD 



1.10573. 



(25) 



So, the indicated inequality is violated. 

At this point, we applied formula (7a) of |27j, giving the Ricci tensor based on the Bures 
metric for diagonal density matrices, 



Ricci(Y, Z) = 3E M)1/) ^- 



V 7 

1 ^L^l ZJ UV 



(26) 



(Pp + P^)(Pm + Pv){p» + Pv) 2 ^(P H + Pv) 
where Y and Z are tangent vectors (traceless Hermitian matrices). We, in fact found, using 
numerical simulations, violations of the lower bound of N for (iV — l)-dimensional manifolds 
on the Ricci tensor required by the Levy-Gromov Inequality, with N = 14 in our case. The 
lowest value we were able to achieve in a series of simulations was 3.45666, so no negative 
values were recorded. (The upper value appeared to be unbounded. We also applied the 



24 



formula (|26|) to the 8- dimensional convex set of 3 x 3 density matrices and found, through 
Monte Carlo simulations, a value of the Ricci curvature as low as 3.00332, so it appears 
conjecturable that 3 is the actual lower bound.) Thus, our evidence here indicates that the 
inequality is not satisfied, apparently since all the conditions for its application have not 
been met. (Using the ansatz elaborated upon earler in sec. II CI we also found numerically 
for the Average and WY monotone metrics that the lower bound of N = 14 on the Ricci 
tensor was violated.) 



1. Area/Volume ratios 

Let us also note here (cf. ly, sec. 6]), in terms of the known values 14j and our con- 
jectures, that for the separable plus nonseparable states, the ratio of the SD 14-dimensional 
hyperarea to the SD 15- dimensional volume, is 

B s + n 8192 

V#= 429^ 6 ' 07831 < 27) 
while for only the separable states, it is 

B ' 8D + PsD = - « 7.53846. (28) 
VSd 13 1 ' 

VIII. CONCLUDING REMARKS 

Needless to say, to the extent any of the conjectures above are, in fact, valid ones, their 
remains the apparently formidable task of finding formal/rigorous proofs. 

A direct /naive "brute force" strategy of symbolically integrating the volume elements of 
the various monotone metrics over the 1 5- dimensional convex sets of separable and all two- 
qubit states — while successful for lower- dimensional scenarios [gj - appears to be quite 
impractical computationally-speaking. It seems that one would have to deal with multiple 
ranges of integration given by fourth-degree polynomials. One might speculate that the 
integration of the (product) measures for the volume elements of monotone metrics over 
the twelve parameterizing Euler angles would yield a result simply proportional to a^ g , 
common to all monotone metrics, and that additional distinguishing factors would appear 
from integrating over the final three variables (8\, 82, #3 in the notation of |4j) parameterizing 
the eigenvalues of the 4x4 density matrices. 
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Perhaps, in this regard, the work of Sommers and Zyczkowski jlj] — which they view "as 
a contribution to the theory of random matrices" — in constructing a general formula for 
V^res f° r -W-level systems, is extendible to (monotone) metrics other than the Bures. The 
volume Vsp n is known, and we have indicated our conjectures that = 64V^ n , V^ g n = 
^&'^Gks = ifso anc ^ Kra?Lr = °°) ^ut we have no similar conjecture, at the present, for 
■ (The Wigner-Yanase metric corresponds to a space of_0.) Bat 
there appears to be no "hint" in the literature as to how one might formally derive simply the 
separable — as opposed to separable plus nonseparable — volumes for any of the monotone 
metrics. 

In this study, we have conjectured that the volumes of separable two-qubit states is, as 
measured in terms of several monotone metrics of interest, simple multiples of the silver 
mean (o~ ag ). It is interesting to point out, it seems, that in certain ( "phyllotactic" ) models 
of the arrangements of (rose) petals, the positions of the petals (in fractions of a full turn) 
are given by the fxactiona. parts of simp.e mu.fp.es of the gM en ratio Q p. 113] Q pp. 
122-123] (cf. |65j, p. 137]). (Of course, it remains possible that we have been somewhat 
"overeager" here to find multiple roles for o^ g . In this regard, a sceptically-inclined reader 
might point out that 2000cr& r ~ 6605.55 is quite close to our sample estimate of 6606.58 
(Table |Vj) for (3-^ Y an d 2ov ~ 6.60555 approximates our sample estimate of 6.60067 for 
B s a - g . Here, a br = is the "bronze mean" Q.) 

In this and other papers , attention has been focused on the matter 

of determining the volumes of quantum states in terms of various monotone metrics. An 
even more considerable body of work concerned with differential geometric properties of the 
monotone metrics is devoted to issues of the scalar curvature of monotone metrics 



5Ja|6J,|67|. 



For instance, it has been found that the scalar curvature of t 
in terms of the Wigner-Yanase metric, | (iV 2 — 1) (iV 2 — 2 



le NxN density matrices Dn is, 



521, while for the Bures metric, it is 



i 1 

24 for N = 2 and for general N, no less than |(5iV 2 — 4)(iV 2 — 1), which value is assumed for 
the fully mixed state 4^1 jf^. (For the m/iniie-dimensional case of thermal squeezed states, 
Twamley 0, eq. (30)] has found the scalar curvature to be given by _ 8 ( cosh ^Z^j^jl ^/ 4 ) ; 
where the "non-unitary" parameter (3 corresponds to the inverse temperature.) It would 
certainly be of interest to find linkages between these two interesting areas of investigation. 
We note that the scalar curvature determines the asymptotic behavior of the volume of a 
Riemannian manifold [3] 0, P- 55, Cor. 5.5 and ex. 3]. Andai 0, eq. (1)] has recently 
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presented a formula for the relation between the volume of a geodesic ball centered at the 
fully mixed state and the scalar curvature there (see also 0, eq. (29)]). 

Our earlier conjecture (j2J) - - in its unadjusted form — as to the exact value of Vg D 
had suggested a similar- type conjecture for qubit-qutrit pairs [12j. Now that we have found 
compelling numerical evidence to reject (and replace it by ()15|)). we obviously must be 
dubious as to the presumed validity of its qubit-qutrit analogue, but presently lack any 
notion as to how to replace it. Additionally, our numerical experience so far indicates 
that it would be extraordinarily difficult to "pinpoint" (accurately estimate) the value of 
the volume of separable qubit-qutrit pairs, since one would then be proceeding in a (more 
computationally demanding) much higher- dimensional (35-d v. 15-6) space, plus the size of 
the separable domain one would be estimating would be much smaller relatively speaking 
(that is, relatively fewer sampled 6x6 density matrices would be separable vis-d-vis the 
4x4 case). 

We summarize in Table IVII our present state of presumed knowledge in regard to the 
various monotone metrics studied here. Of course, one would aspire to find the junctionals 
that map an operator monotone function f metr ic(t) into V^ etric , V^ t n ric , B s metHc1 B s +£ ric and 

metric- 
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TABLE VI: Conjectured values (except for V* + - n and B s +™ , which are known) oiV s - . , V s+ ~ n , 

v Bures Bures ' metric metric 

£s+n rs an( j ft ~ for (four times) various monotone metrics, listed in order of increasing 

metric' metric ^metric V ) J to 

volume size, together with the corresponding operator monotone functions and Morozova-Chentsov 
functions. (For the various denominators, we have the interesting prime decompositions: 5040 = 



2 4 • 3 2 • 5 • 


7; 1750 = 


2 


• 5 3 • 7; 8448 


= 2 8 • 3 


• 11; 315 = 


= 3 2 • 5 • 


7; 135135 


= 3 3 • 5 • 


7 • 11 • 13; 


42075 = 3 


2 • 5 2 • 11 • 


17 and 45045 = 


= 3 2 • 5 • 


v.ii. 


13. 


As pertains to numerators: 


512 = 2 9 ; 


7735 = 5 • 


7 • 13 • 17; 


262144 = 2 18 ; 


15950 = 


2 • 5 2 


• 11 


• 29; 255 


= 3 ■ 5 ■ 


17; 495 = 


3 3 • 5 • 11; 


270 = 2 • 3 


3 • 5; and 616 


= 2 3 -7 - 11.) 
















metric 


/(*) 






V s 


ys+n 




B s 




P 


B s + f3 


Bures 


l+t 




2 


<JAg 


7T 8 




A3aA g 


5127T 7 


55<TAg 


98a Ag 


2 




Pp+py 


3 


5040 




39 


135135 


39 


39 


GKS 


t t/(t-l) 
e 




Pn 

e (Pj±) Pv-Pp 
■ Pv ' 

Pv 


5 


1750 




? 


? 


270a A9 
77 


? 


WY 


(Vi+i) 2 




A 




? 


7735cta 3 - 


262144tt 7 


l5950a A g 


23685(7A 5 


A 




(Vpv+Vpv) 2 


4 




45045 


Avg 


l+6t+t 2 




Mpp+Pv) 




257T 8 




255cta 9 


3437tt 7 


lh(JAg 


495cta 9 


A+At 




p'l+§pp.pv+p'i 


9 


8448 




16 


42075 


16 


KM 


t-1 




log^ 

b Pv 


Wa Ag 


47T 8 




OO 


oo 


616(ta 9 


oo 


log* 




Pp-Pv 


315 




13 


NI 


2(t-l) 2 




(p M +p.)log 2 ^ 








OO 


oo 


7 


oo 


(l + t)(logt) 2 














maximal 


2t 






oo 


OO 




OO 


oo 


oo 


oo 



[5] K. Zyczkowski, P. Horodecki, A. Sanpera, and M. Lewenstein, Phys. Rev. A 58, 883 (1998). 
[6] K. Zyczkowski, Phys. Rev. A 60, 3496 (1999). 
[7] P. B. Slater, J. Phys. A 32, 5261 (1999). 
[8] P. B. Slater, Euro. Phys. J. B 17, 471 (2000). 
[9] R. F. Werner, Phys. Rev. A 40, 4277 (1989). 
[10] A. Peres, Phys. Rev. Lett. 77, 1413 (1996). 

[11] M. Horodecki, P. Horodecki, and R. Horodecki, Phys. Lett. A 223, 1 (1996). 
[12] P. B. Slater, J. Opt. B 5, S691 (2003). 

[13] J. F. Verstraete and B. D. Moor, Phys. Rev. A 64, 010101 (2001). 
[14] H.-J. Sommers and K. Zyczkowski, J. Phys. A 36, 10883 (2003). 
[15] P. B. Slater, J. Phys. A 32, 8231 (1999). 

[16] K. Zyczkowski and H.-J. Sommers, J. Phys. A 36, 10115 (2003). 

[17] T. Tilma, M. Byrd, and E. C. G. Sudarshan, J. Phys. A 35, 10445 (2002). 



28 



[18] J. N. Lyness, Math. Comput. 19, 394 (1965). 
[19] G. Okten, MATHEMATICA in Educ. Res. 8, 52 (1999). 
[20] D. Petz and C. Sudar, J. Math. Phys. 37, 2662 (1996). 
[21] D. Petz, Lin. Alg. Applies. 244, 81 (1996). 

[22] A. Lesniewski and M. B. Ruskai, J. Math. Phys. 40, 5702 (1999). 

[23] R. E. Kass, Geometrical Foundations of Asymptotic Inference (John Wiley, New York, 1997). 
[24] M. Hiibner, Phys. Lett. A 63, 239 (1992). 
[25] M. Hiibner, Phys. Lett. A 179, 226 (179). 
[26] J. Dittmann, J. Phys. A 32, 2663 (1999). 
[27] J. Dittmann, J. Geom. Phys. 31, 16 (1999). 

[28] H. P. Yuen and M. Lax, IEEE Trans. Inform. Th. 19, 740 (1973). 
[29] H. Hasegawa, Rep. Math. Phys. 39, 49 (1997). 
[30] D. Petz, J. Math, Phys. 35, 780 (1994). 

[31] P. W. Michor, D. Petz, and A. Andai, Infin. Dimens. Anal. Quantum Probab. Relat. Top. 3, 
199 (2002). 

[32] M. R. Grasselli and R. F. Streater, Infin. Dimens. Anal. Quantum Probab. Relat. Top. 4, 173 
(2001). 

[33] H. Faure and S. Tezuka, in Monte Carlo and Quasi-Monte Carlo Methods 2000 (Hong Kong), 

edited by K. T. Tang, F. J. Hickernell, and H. Niederreiter (Springer, Berlin, 2002), p. 242. 
[34] G. Call and D. Velleman, Amer. Math. Monthly 100, 372 (1993). 
[35] H. S. Hong and F. J. Hickernell, ACM Trans. Math. Software 29, 95 (2003). 
[36] J. Matousek, Geometric Discrepancy: An Illustrated Guide (Springer, Berlin, 1999). 
[37] M. J. W. Hall, Phys. Lett. A 242, 123 (1998). 
[38] G. A. Christos and T. Gherghetta, Phys. Rev. A 44, 898 (1991). 
[39] S. Abe and A. K. Rajagopal, Phys. Rev. A 60, 3461 (1999). 

[40] R. D. Gill, G. Weihs, A. Zeilinger, and M. Zukowski, Proc. Natl. Acad. Sci. 99, 14632 (2002). 
[41] V. W. de Spinadel, Int. Math. J. 2, 279 (2002). 
[42] G. Gumbs, J. Phys. A 22, 951 (1989). 

[43] J. Kappraff, Beyond Measure (World Scientific, River Edge, NJ, 2002). 
[44] J. G. Escudero and J. G. Garcia, J. Phys. Soc. Japan. 70, 3511 (2001). 
[45] M. S. El Naschie, Chaos, Sol. Fract. 10, 1303 (1999). 

29 



[46] W.-F. Chuan and F. Yu, Fibonacci Quart. 38, 425 (2000). 

[47] C. Bavard and P. Pansu, Ann. Sci. Ecole Norm. Sup. 19, 479 (1986). 

[48] B. H. Bowditch, J. Austral. Math. Soc. Scr. A 55, 23 (1993). 

[49] R. P. Bambah, V. C. Dumir, and R. J. Hans-Gill, Studia Sci. Math. Hungar. 21, 135 (1986). 

[50] C. Krattenthaler and P. B. Slater, IEEE. Trans. Inform. Th. 46, 801 (2000). 

[51] P. B. Slater, J. Phys. A 34, 7029 (2001). 

[52] P. Gibilisco and T. Isola, J. Math. Phys. 44, 3752 (2003). 

[53] E. Wigner and M. Yanase, Proc. Natl. Acad. Sci. 49, 910 (1963). 

[54] S. Luo, Phys. Rev. Lett. 91, 180403 (2003). 

[55] P. B. Slater, Phys. Lett. A 247, 1 (1998). 

[56] S. L. Braunstein, C. M. Caves, R. Jozsa, N. Linden, S. Popescu, and R. Schack, Phys. Rev. 

Lett. 83, 1054 (1999). 

[57] R. Clifton and H. Halvorson, Phys. Rev. A 61, 0121018 (2000). 

[58] L. Gurvits and H. Barnum, Phys. Rev. A 68, 042310 (2003). 

[59] S. Szarek, quant-ph/0310061. 

[60] F. J. Bloore, J. Phys. A. 9, 2059 (1976). 

[61] M. Shi and J. Du, quant-ph/0103016. 

[62] M. Gromov, Metric Structures for Riemannian and Non-Riemannian Spaces (Birkhauser, 
Boston, 1999). 

[63] M. Livio, The Golden Ratio (Broadway, New York, 2002). 

[64] J. H. Conway and R. K. Guy, The Book of Numbers (Springer- Verlag, New York, 1996). 

[65] M. H. Freedman, Comm. Math. Phys. 234, 129 (2003). 

[66] A. Andai, J. Math. Phys. 44, 3676 (2003). 

[67] J. Dittmann, J. Geom. Phys. 31, 16 (1999). 

[68] J. Twamley, J. Phys. A 29, 3723 (1996). 

[69] T. Sakai, Riemannian Geometry (Amer. Math. Soc, Providence, 1996). 

[70] A. Andai, math-ph/0310064. 

[71] D. Petz, J. Phys. A 35, 929 (2002). 



30 



